Normal Mode Spectrum of the Deuteron in 
the Skyrme Model 



C. Barnes^, W. K. Baskerville^'^ and N. Turok*^ 



. . . '^Department of Physics, Princeton University, Princeton, NJ 08544-0708 USA 

Q>^ ■ ^Physics Department, University of Wales Swansea, Singleton Park, 

G\ '. Swansea SA2 8PP 



'^Department of Applied Mathematics and Theoretical Physics, Silver Street, 
Dh! Cambridge CB3 9EW UK 

< 

^ ■ Abstract 
> 
OO 

' The normal mode spectrum of the deuteron in the Skyrme model is computed. We 

■ find a bound doublet mode below the pion mass, which can be related to the well- 

known 90° scattering of two skyrmions. We also find a singlet 'breather' mode and 
I another doublet above the pion mass. The qualitative pattern of the spectrum is 

On ' similar to that recently found for the B = A multiskyrmion. The symmetries of all 

the vibrational modes are presented. 



^ ■ 1 Introduction 



The Skyrme model [1] has had some success in describing both single nucleon 
properties and the nucleon-nucleon interaction [2-4] . It is therefore of interest 
to see how well the model performs for larger nuclei. Classical multiskyrmion 
solutions are known up to baryon number 5 = 9 [5,6]. However, these must 
be quantised before any comparison to the real world can be made. 

It is generally agreed that a proper treatment of the Skyrme model as a quan- 
tum field theory is difficult. Instead, a limited number of degrees of freedom 
are quantised, and usually only the zero mode collective coordinates [2]. This 
reduces the model to finite-dimensional quantum mechanics. The conventional 
wisdom is that 6A^ degrees of freedom are required to describe a system con- 
taining nucleons (as would be required for widely separated skyrmions). 
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However, minimal energy solutions for 5 > 1 are in general single large soli- 
tons, which have a maximum of 9 zero modes. Therefore even for a very limited 
quantisation it is necessary to include additional discrete modes. 

This is the principal motivation for the calculation of the vibrational spectra of 
multiskyrmions. A full quantisation should also include the effects of the soli- 
ton on the vacuum fluctuations, or radiation. We shall not attempt this here, 
but will merely present a new method for computing normal mode spectra, 
with particular reference to the B — 2 solution. The results ior B — A are de- 
scribed in a concurrent paper [9] . A similar calculation was recently performed 
by Walet [7] for 5 = 2 and B = 3, using the Yang-Mills instanton-induced 
form for the Skyrme fields. The main differences in the present calculation 
is that we do not employ this approximation, and also explicitly include a 
pion mass term in the Lagrangian. This mass term clearly plays a key role in 
determining which vibrational modes are bound. 

The vibrational spectra are interesting in their own right, especially since 
we have discovered hitherto unsuspected similarities between the two cases 
considered. In particular, the spectrum divides into two types of modes. The 
lower modes correspond to those expected from an approximate correspon- 
dence between BPS monopoles and Skyrmions [9]. Going up in frequency, one 
encounters first the 'breather' and then higher multipole breathing modes. 
While the structure of the spectrum is interesting, we emphasise that it repre- 
sents merely a first step in a rather more ambitious long term project, namely 
the computation of properties of quantised multiSkyrmions. 



2 Computation of the Spectrum 

The basic idea behind our method for calculating the bound state spectrum is 
quite simple. First, the classical solution is obtained by numerical relaxation. 
It is then very slightly perturbed, and the resulting configuration evolved for a 
long time at constant energy. During this process, field values at certain points 
are recorded as a function of time. This data can then be Fourier analysed to 
obtain power spectra. Both the relaxation to the classical solution, and the 
subsequent time evolution after perturbation, require numerical solution of 
the Skyrme equations. The development of a fast computer code to perform 
this task lies at the heart of our method. In the current letter, we shall only 
briefly outline the main ideas underlying our algorithm; further details may 
be found in [10]. 

The Skyrme model has three free parameters: Fj^, e and m^r. The work of 
Walet [7] indicates that binding energies after quantisation probably depend 
very sensitively on these parameters. We defer further discussion of parameters 



2 



until the ftill calculation has been performed. Instead, we will express our 
results in terms of the usual dimensionless Skyrme units [12]. In these units, the 
topological lower bound on the energy is 127r^B, and the Skyrme Lagrangian 
density is 

jC = lTr{L,L^) + ^Tri[L^,L,][L^,L''])+Ap\TrU - 2), (1) 

where /3 = rriT^/FT^e is a dimensionless constant, = 17^8^17 and U = a + m.f 
is the SU(2)-valued scalar field (tj are the Pauh matrices). For definiteness, 
we have set (3 — 0.263, following Adkins and Nappi[8], although eventually we 
would would hope to extend the calculation to treat /3 as a free parameter. 

The action S = J d'^C is discretised on a finite lattice, and the equations of 
motion are derived by varying with respect to each field component j, k, t) 
(where (j) = (cr, ff)), at each point in the lattice. The SU{2) constraint is 
enforced by the introduction of a term \{(f)'^ — 1), where A is a Lagrange 
multiplier. At this point we run into the usual problem that the "kinetic" (ie. 
time dependent) part of the quartic Skyrme term causes coupling between 
nearest neighbours in time and space. The problem is trivially solved for the 
initial relaxation by simply ignoring the offending terms, which are irrelevant 
for a static solution. But to compute the normal modes we must consider small 
fluctuations about the static solution: = 0o + S(f). If we linearise in 6(f) (just 
for the problematic terms), then the system of equations uncouples, and we 
are left with an equation of the form 

j, k,t+l)^~Xr{^, J, k, t) + R%i, J, k, t) (2) 

at each point, where A is a constant related to the Lagrange multiplier. The 
SU{2) constraint can then be invoked for 0(t+l), to produce a quadratic 
equation which can be solved for A. This can then be substituted back into 
Equation (2) to define the time evolution. 



3 Results and Interpretation of Spectra. 

Our simulations are carried out on a finite lattice of A^'^ points, of total volume 
L^, so that the lattice spacing Ax — L/N. Deuteron spectra have been calcu- 
lated for a variety of box sizes and latice spacings. We find that the frequencies 
observed show no discernible dependence on lattice spacing, provided that the 
lattice spacing is no larger than about 0.125 in our length units. There is, 
however, significant dependence on box size. Figure 1 shows deuteron spec- 
tra for L = 6 and L — 8. Our code uses periodic boundary conditions, so 
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Fig. 1. Deuteron power spectra, for two different box sizes (L=6 and L=8). 
Log(power) is plotted against angular frequency uj. 

the finite size effects can be understood in terms of the interaction between 
neighbouring 'image' solitons. Since a pion mass term is included explicitly in 
the Lagrangian, these effects should decay exponentially as L increases. The 
frequencies observed for L = 8 should already be a reasonable approximation 
to the infinite separation (box size) limit; however, this can eventually be im- 
proved, with more data, by fitting the exponential curve, and extrapolating. 

Having obtained the normal mode spectrum, the next problem is to interpret 
it. It is necessary to distinguish between vibrational modes and radiation, and 
we also wish to determine the degeneracy and symmetry of all modes. Note 
that radiative modes are discrete, rather than forming a continuum above the 
pion mass, because of the finite box size. For the same reason, bound states are 
possible above the pion mass. In an infinite volume box, radiation continuum 
modes at the frequencies of the vibrational modes would mix with the latter, 
causing them to be spatially extended out to infinity. In a finite box, however, 
the relevant continuum modes may be absent, causing the vibrational modes 
above the pion mass to be spatially localised. 

To gain information about the observed peaks, the process of time evolution 
is repeated, this time projecting out all previously identified frequencies. We 
thus obtain Fourier amplitudes for each mode, at all points on the lattice for 
each field component 

Al{i,j,k) = ydt50'^(i,j,k,t)cos^t (3) 



where 50" = 0" — 0g is the difference between the actual field value at a given 
time, and the static solution. For each frequency, the total time of integration 
is chosen to be an exact multiple of the period of oscillation. 
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The inner product between two such amphtudes can be defined 

{A^\A^,) = j dx^ K^'^(x) A:(x) A^,(x), (4) 

with summation imphed, where K"^ — (5"''(l + (9i0o)^) — ^i^o'^i^o matrix 
which multiphes d^cf)'^ on the left hand side of the (hneariscd) Skyrme equa- 
tions. It effectively plays the role of a spatially dependent inertia tensor in the 
dynamics. The inner product is defined in this way to ensure orthogonality 
between modes of different frequency. 

The most immediate concern is to distinguish between bound modes and ra- 
diation in the spectrum. To do this, it is useful to calculate the norm of a 
given mode as a function of radius. Instead of integrating {A^\A^) over the 
whole box, we instead integrate only for points in a cube in the middle. By 
varying N, {A^\AJ) is mapped out as a function of "radius" (or side length). 
If a particular mode is bound, then it should be localised at or near the static 
solution. Thus the curve of (A^^IA;^) vs should rise steeply at a radius cor- 
responding roughly to the soliton, then fiatten out as the edge of the box is 
approached. Radiation, on the other hand, should be spread evenly through- 
out the box, giving a curve going approximately as N^, rising most steeply 
at the edge. The curves for the six modes labelled in Figure 1 are plotted in 
Figure 2. It is clear that the curves do indeed display characteristic shapes, 
so that this method successfully provides a clear distinction between bound 
modes and radiation. We conclude that modes 3 and 4 are radiative, whereas 
modes 1, 2, 5 and 6 are bound. 

The next obvious question is the degeneracy of all the modes. If a given fre- 
quency u) is n-fold degenerate, then an mxm matrix of inner products (74^|74^) 
(where m > n), will only have n non-zero eigenvalues. Fourier amplitudes from 
different runs, with different random initial perturbations, (but for the same 
frequency), can be used to calculate the degeneracies via this method. How- 
ever, cleaner results are often obtained by taking one amplitude, then gener- 
ating more by applying different symmetry operations. The symmetries used 
must leave the static solution invariant; they generally consist of a physical 
transformation followed by an isospin transformation of the pion fields. Using 
this method, we find that modes 3 and 5 are singlets, while peaks 1, 2, 4, and 
6 correspond to doublets. 

The only remaining issue concerns symmetry. We wish to classify all the modes 
according to the representation they form of the symmetry group of the static 
solution. The classical B = 2 multiskyrmion is axially symmetric, and in addi- 
tion possesses a refiection symmetry in the plane of the torus. The symmetry 
group is -Doo/i; axial symmetry extended by inversion. Its character table is 
given in Table 1. The notation of Hamermesh[ll] has been used for the conju- 
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Norm of Mode vs Radius: L=6, N=49 
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Fig. 2. VS N, for the six lowest modes in the deuteron spectrum. All modes 

are normalised so that {A^\AJ) = 1 when integrated over the whole box 

gacy classes and the singlet representations. E is the identity, C(0) is rotation 
by angle (p, is reflection in any plane including the z-axis, and / is inversion. 
The singlet representations are labelled or to indicate whether they are 
even or odd under parity, while the + or — refers to their behaviour under a^. 
There are an infinite number of doublets, labelled by an integer n, and parity. 



It is clear that singlet representations are distinguished by their properties 
under (y.^, and I. The symmetries of modes 3 and 5 can therefore be determined 
by taking matrix elements of these operators 
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Table 1 

Character table for £)ooh 
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(3|a,|3) = 1 (3|/|3) = -1 

(5|(7„|5) = 1 (5|/|5) = 1. (5) 

For doublet states, the simplest way to determine n is to search for the lowest 
angle d such that (X|C(0)|X) = (this implies that (X|C(20)|X) = -1 and 
{X\C{A(j))\X) = 1). If the lowest angle = 90° then n = 1, if = 45° then 
n — 2, and so on. In practice, n = 1 or 2 are the only realistic values, so 
the matrix elements of C(45°) and C(90°), together with inversion, provide 
sufficient information to distinguish the representations. 



(1|C(45°)|1) 


1 


(1|C(90°)|1) = 


(1|J|1) = 1 


(2|C(45°)|2) 


= 


(2|C(90°)|2) = -1 


mm = 1 


(4|C(45°)|4) 


= 


(4|C(90°)|4) = -1 


(4|7|4) = 1 


(6|C(45°)|6) 


1 

~ V2 


(6|C(90°)|6) = 


(6|7|6) = -1 



To summarise, the first peak is a doublet, with symmetry 1+. This is the 
representation under which axial vectors transform in D^h- These states do 
not correspond to true vibrations, but rather represent two broken zero modes, 
spatial rotations of the torus around the x and y axes. The periodic boundary 
conditions of the lattice introduce the possibility of low energy 'spin waves', 
whose frequency will decrease as the box size increases, eventually reverting to 
zero modes in the infinite limit. This explains the extreme shift in frequency 
observed for this mode when the box size is changed from L — 6 to L — 8. 

The second peak in the spectrum represents the first true finite frequency 
bound modes. They are doubly degenerate, with a quadrupole symmetry 2+. 
If amplified, these would become the classic 90° scattering of two skyrmions. 
There is also a close analogy to the scattering of BPS monopoles. 

The zero modes of the two-monopole toroidal BPS solution may be fairly 
simply understood. The BPS solutions are associated with rational maps con- 
structed as follows. A complex variable z is associated with each direction in 
space through a projective map. The complex z plane is embedded in three 
Euclidean dimensions, and a unit two-sphere is drawn centered on the origin. 
Each direction n defines a point on the two-sphere, which is then projected 
from one of the poles of the two-sphere onto the z plane. One then defines 
a rational map R which is a ratio of polynomials in z, each up to z'^ for a 
N-monopole solution. Just as z defines a direction in real space, R defines 
one in internal space, acted on by the global SU (2) symmetry via Moebius 
transformations R — >• {aR + /?)/(— -|- a*), with a and (3 being the 11 and 
12 matrix elements of the SU{2) matrix. 
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For special choices of i?, the solutions exhibit enlarged symmetry. In our 
case, the 2-monopole solution represented hj R = is the toriodal solu- 
tion which has -Doo/i symmetry. It is invariant under rotations z — > e^'^z com- 
bined with isospin rotations R — > e~^"^it!, and also inversions z — > —l/z*, 
R — > 1/R*. The general perturbed 2-monopole solution is represented by 
R = {z'^ + az + b) I [cz^ + (i^; + 1 -|- e) where the coefficient of z^ in the numer- 
ator is taken to be unity and the complex parameters a, 6, c, (i, e are small. 
These parameters are not invariant under D^^^ but transform according to 
some irreducible representation, which has ten real dimensions. It is straight- 
forward to check by computing the characters of the symmetry group elements 
that the representation is 2+ -|- 2~-|- 1+ -|- S~ -|- 1~ -|- E~. The isospin zero 
modes comprise the 2~ representation, the rotational zero modes comprise 
the 1+ + Eg representation and the tranlational zero modes the 1~ -|- S~ 
representation. This leaves the 2+ representation as the remaining nontrivial 
zero mode of the 2-monopole BPS solution. And this representation coincides 
with our results for the lowest vibrational mode of the deuteron in the Skyrme 
model. 

The next two frequencies correspond to radiation, with wave number k — 
0. These modes are spht, rather than forming a degenerate triplet, because 
the pion fields of the static solution do not all transform the same way. The 
radiative modes fall into the same representations as the static solution pions: 
+ 2+. 

The fifth peak is a breathing mode, where the size of the soliton fluctuates 
with no change in shape. It is a singlet, and transforms trivially under -Doo/i, 
as would be expected. The sixth peak represents the last bound states. These 
form a doublet, with symmetry 1~. The torus expands on one side while 
contracting on the other, a dipole 'breathing' motion. Altogether then, we 
have found five finite energy normal modes in the deuteron spectrum. 



4 Conclusions. 

The normal mode spectrum of the deuteron has been calculated. The nature 

of all states has been identified, and their degeneracies and symmetries calcu- 
lated. We find a bound doublet below the pion mass, which is related to 90° 
scattering between two skyrmions. Two more bound states appear above the 
pion mass: a breather, plus a second doublet with a dipole symmetry. 

The overall pattern of the deuteron spectrum displays a remarkable similarity 
to that recently found for the 5 = 4 solution [9]. There again, all modes 
below the pion mass can be related to those expected from an approximate 
correspondence with the zero modes of BPS multimonopoles. The first B — ^ 
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bound mode above the pion mass is a breather, and the remaining bound 
modes involve 'breathing' motion of the baryon density. It is tempting to 
speculate that a similar pattern will be observed for all multiskyrmion spectra, 
but more examples are required to confirm this. 

One further remark may be worthwhile. We have found five finite energy bound 
modes in the deuteron spectrum. Added to the eight zero modes, this makes 
a total of 13. Sixteen bound modes were found for B = 4, giving a total of 25. 
In both cases, then, there is a total of 6B + 1 modes, one more than is usually 
expected. The "extra" mode is the breather; every multiskyrmion, including 
the single skyrmion, must have such a mode. The remaining vibrational modes 
can be considered as corresponding to the broken zero modes of the separated 
skyrmions. 
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Fig. 3. Deuteron vibrational modes: plots of B{(j)o + S(l)) — B{(pQ). Isosurfaces indicate 
surfaces of constant negative perturbation; contours positive. Note that all modes 
except mode 3 are plotted looking more or less directly down the z-axis, so that 
whereas modes 2 and 4 are quadrupoles in the (x, y)-plane, mode 1 is a dipole, but 
is odd under reflections in this plane. Mode 3 is odd under z-reflections, but the 
twisting is merely an artifact due to slight contamination from mode 2. 
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